How to Measure the Quantum State of Collective Atomic Spin Excitation 
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The spin state of an atomic ensemble can be viewed as two bosonic modes, i.e., a quantum 
signal mode and a c- numbered "local oscillator" mode when large numbers of spin- 1/2 atoms are 
spin-polarized along a certain axis and collectively manipulated within the vicinity of the axis. We 
present a concrete procedure which determines the spin-excitation- number distribution, i.e., the 
diagonal elements of the density matrix in the Dicke basis for the collective spin state. By seeing 
the collective spin state eis a statistical mixture of the inherently-entangled Dicke states, the physical 
picture of its multi-particle entanglement is made clear. 

PACS numbers: 03.67.Mn, 42.50.Dv, 42.50.Fx 



I. INTRODUCTION 

Squeezing the spin projection noise 1] below the stan- 
dard quantum limit [2] necessarily means multi-particle 
entanglement among elementary spins 0, Q. Such a 
squeezed spin state (SSS) H is believed to improve the 
precision of atomic clocks UIE Qjand to boost the sen- 
sitivity of atomic magnetometers and atomic interfer- 
ometers [ollioj. 

The SSSs have been generated in several laborato- 
ries since the pioneering experiment with cold cesium 
atoms [Tl| . To characterize the quantum nature of the 
SSSs, the spin operators are usually bosonized such that 
Xa = Jy/VJ^ and Pa = Jz/VJ^-, where [Jy, Jz] = iJx 
with Jx = ( Ja;) being c-numbered quantity [Til |12J and 
these variances are analyzed in an analogous fashion 
to quantum optics Given that in quantum optics 

single-mode squeezed state does not evoke the concept of 
multi-particle entanglement, the connection between the 
SSSs and their entanglement seems to be obscure at first 
glance. 

Here we provide a different view on the collective spin 
state. In our framework the spin operators are considered 
to be made up of two bosonic modes, i.e., a quantum sig- 
nal mode and a c-numbered "local oscillator (LO)" mode. 
This viewpoint smoothly links the quadratures and the 
density matrix in the inherently-entangled Dicke basis, 
just like the link between the quadratures and the den- 
sity matrix in the inherently-nonclassical Fock basis in 
quantum optics 0, 0, 0| . We show a concrete proce- 
dure to reconstruct the spin-excitation-number distribu- 
tion, i.e., the diagonal elements of the density matrix in 
the Dicke basis for general collective spin states including 
the SSSs. The resultant statistics of spin excitations give 



us clear physical picture of multi-particle entanglement 
for collective spin states. 



II. SYMMETRIC COLLECTIVE SPIN AS TWO 
BOSONIC MODES 

We start by examining the optical homodyne detection 
from a perspective more appropriate for grasping the con- 
cept of "LO" mode in a collective spin system. Let by 
and bn denote the annihilation operators for the quan- 
tum mode and the LO mode, respectively. The dynam- 
ics of the two bosonic modes can be analyzed by angular 
momentum operators (Stokes operators) [T3 |. which are 
given by 



Sx — 

Sy — 

= 



{Vybv - Vh^h) 

{bVhH + h^Hbv) 
-i{VybH - b^H^\ 



(1) 
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The operator bn for the LO mode can be assumed 
to be a c-numbered quantity ^fne~^'^ since the num- 
ber of photons excited in this mode, {bl^bn) — n, is 
enormously large. The angular momentum operator Sz 
then becomes Sz = —V2n(smip qv + cos Pv), where 
qv = {pv + ^v)/V2 and pv = iiby — 6y)/\/2 are the 
quadrature operators. The operator Sz is thus nothing 
but the observable of the homodyne detection with 
ip representing the LO phase and y/n representing the 
mean LO amplitude. 

The converse of the above argument also holds: the 
dynamics of a collection of spin-1/2 atoms with the 
symmetric collective spin operators Jx^ Jy, and Jz, 
can be analyzed by two bosonic modes. The correspond- 
ing annihilation operators and 0| for spin-up and spin- 
down modes are defined as a||M)| = ^/M\AI — l)t and 
ai\M)i = ^/M\M — respectively. Since the x-axis is 
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defined as the quantization axis, we have 



J:r. = 



J 11 —- 



1 



At; 



J. 



(2) 



When all atoms are condensed in the spin-down mode, 
the state can be written as |0) | (g) | A^) | . This state is called 
the coherent spin state (CSS) .2] which defines the stan- 
dard quantum limit, and is the eigenstate of the operator 
in Eq. © with the eigenvalue -N/2. |M)| (g) \N-M) i 

then corresponds to the eigenstate of the operator 
with the eigenvalue —N/2 + M, that is, M atoms excited 
into the spin-up mode with N — M atoms remaining spin- 
down, and is a multi-particle entangled state except for 
M = and N. For instance, the state (g) |iV - 1) i 
can be described in terms of elementary spins (| 1)i and 
I as 



|l)T®l^-l)i 



-^(IT>iU>.^ 



-^(U>i|T>.^ 



-7^(U>iU>2 



(3) 



and does have multi-particle entanglement 4J. Given 
that a collective state belongs to the symmetric subspace 
of the full Hilbert space (spanned by 2^ states of N spin- 
1/2 atoms), it can be expressed by using these entangled 
states, namely, by the Dicke basis. 

With this second quantized formalism of spin degree 
of freedom, the above treatment of optical homodyne de- 
tection can be translated into the language of a collective 
spin. For states in the vicinity of the CSS (|0)| |iV)|), 
the operator can be modified to 



Jz 



- v/iV/2(sin d + cos p| ) 



(4) 



This follows from the assumption that the operator for 
the spin-down mode is a c- numbered quantity ViVe"*'', 
since the number of spin-down atoms ((ajaj) — N) is 
enormously large. Here g| = (6|-|-&|)/-\/2 andp-f = i(6| — 

6|)/-\/2 are the collective spin analogues of the quadra- 
ture operators (the physical meaning of the quadrature 
distribution for the collective spin is briefly discussed in 
the Appendix). This approximation is similar to that 
for the optical homodyne detection and somewhat dif- 
ferent from the spin quadratures Xa and Pa 0, 0| as 
we explicitly introduce two bosonic modes. The atoms 
condensed in the spin-down mode act as the "LO" mode, 
collectively amplifying the tiny differences of the quadra- 
ture distributions for the spin-up atomic states. This sug- 
gests that tomographic measurements [T^ of the quan- 
tum state of the collective spin excitations, such as the 



squeezed spin state (SSS) "s"] and the Dicke state [l^ . 
may be possible in an analogous fashion to the tomog- 
raphy of the squeezed state of light 0| and the single 
photon state respectively. 



III. TOMOGRAPHIC MEASUREMENT OF 
COLLECTIVE SPIN 

Let us consider how to measure the quadrature for a 
collective spin, i.e., the operator Jz in Eq. 10} . In order 
to do this, it is required us to access only the symmetric 
collective spin operators i.e., to vahdate the second 
quantized formalism of spin. Although the most intuitive 
approach may be to measure Jz by the Stern-Gerlach 
technique, the atoms respond to Stern-Gerlach magnetic 
field individually rather than collectively 18j. Kuzmich 
et al. showed that the atomic ensemble collectively inter- 
acts with an off-resonent optical probe and that the col- 
lective spin operator Jz can be measured in a quantum 
nondemolition (QND) way Recently, Julsgaard et 

al. ^3 have applied this technique to the tomography 
of collective spin in the context of the quantum mem- 
ory experiment, in which they proved that the conjugate 
variables of light, i.e., the quadratures, and Pl, were 
faithfully mapped onto the atomic quadratures Pa and 
Xa (pt and in our language), respectively. The ef- 
fective Hamiltonian describing spin QND measurement, 
namely dispersive vector light-shift interaction, is given 
by J^,21,] 



Hi = hgSzJz, 



(5) 



where g is the field-atom coupling coefRcient, Sz is the 
Stokes operator defined by Eq. and Jz is the angular 
momentum operator defined by Eq. (|2Jl. Note that the 
bosonic operators by and bn in Eq. now represent the 
annihilation operators for horizontally and vertically po- 
larized photons, respectively. In the Heisenberg picture, 
the set of Stokes operators, which is initially defined by 
S(0) = {5'^(0), Sy{0), Sz{0)}, evolves into 







Sy{r) 




I Sz{r) \ 





cosigrJz) -sinigrJz) 
ain{gTJz) cos(gTj^) 
1 



Sy{0) 

Sz{0) 

(6) 

after passage of the interaction time r. The set of an- 
gular momentum operators, J(0), evolves similarly with 
J and S interchanged. Here, Jz(0) = Jz{t) = Jz and 
5^(0) = Sz{t) — Sz since Jz and Sz commutes with 
the interaction Hamiltonian given by Eq. ||SJ). When 
the initial state for the optical probe is set for horizon- 
tally polarized light, |0)y (g) \n)H, which is given by the 
eigenstate of the operator Sx{0) with the eigenvalue —n. 
Then the measurement of the operator S ,,(t) , which can 
be realized by the balanced polarimeter |l9j | , effectively 
works as the quadrature measurement for a collective 
spin. To illustrate this concept, the operator Sy{T) in 
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Eq. ® 

can be rewritten as S'y(r) « gTJz{0)Sx{0)-\- Sy{0) 
by assuming gT{Jz) <C 1. Considering the initial con- 
dition for the optical probe, the operators Sx{0) and 
Sy{0) can be approximated by Sx{0) — —n and 5j^(0) = 
■\/n(6y(0) -I- &\/(0)) = \/2rigy(0), respectively. This ma- 
nipulation again applies the fact that the operator 6ij(0) 
is deemed to be a c-numbered quantity ■y/ne"*'^, that is, 
(6|^(0)6_f/(0)) = n » 1. Here, the phase Lp between the 
two modes, |0)y and In)//, is irrelevant and is set at 0. 
We then have 

(t) w gTn^/Nj2{sm g| (0) + cos p| (0) ) -f \^iqv (0) . 

(7) 

The "LO" phase for the collective spin, corresponds 
to the azimuthal angle about the x-axis and can be con- 
trolled by the magnetic field along the x-axis 0, [23, 
|2^ . The observable of the balanced polarimeter, Sy{T), 
can thus be viewed as the amplified quadrature for a 
collective spin, though it includes an additional term 
'\/2ngy(0), which represents the shot noise due to the 
optical probe. 

To keep things simple let us normalize the operator 
Sy{T) as Q EE Sy{T)/M where M = {{grn^/Nplf + 

{V2nf) ' . Then we have 

Q - V^(sin79 gt(0) +C0SZ9 pt(0))V^^9v(0) (8) 

with ^ = gTn^/Nj2/J\f. Let P(Q,i?) = 
(Q, z9|(|0)yy(0|(8)p|)|(3, 1?) denote the histogram obtained 
by measurement of the operator Q, where IQ,")?) is the 
eigenstate of the operator Q with the eigenvalue Q. 
Here, p| represents the density matrix of the spin-up 
mode before field-atom interaction. The optical-noise 
term •\/2ngv(0) prevents simple reconstruction of the 
density matrix p-^. This is essentially the same situa- 
tion where detector inefficiencies and optical losses de- 
grade the quadrature distribution for the original pho- 
tonic state, necessitating loss-error compensation for re- 
construction of the original state [TsL 123 . . Following 
the argument on detection-loss induced error in the op- 
tical homodyne detection the histogram V[Q,d) is 
given by 

P(Q,^) = r dq' P;^)(<z') Pv(^^^), (9) 

with the quadrature distribution for the collective spin 
state 

¥\'>\q)^^{q,d\p^\q,{)h, (10) 

and the i?-invariant quadrature distribution for the opti- 
cal probe 

Pv(g)=y(q|(|0)yv(0|)|q)y = (l/V^)exp(V),, (H) 

which is the Gaussian distribution discribing the shot 
noise. Here |(7, i?)t and \q)v are the eigenstates of the the 



operators (sini? q^ +cos'd p^) and qv^ respectively. Mea- 
surement of Syij) thus yields the smoothed quadrature 
distribution for the collective spin. The tomographic re- 
construction of /3| can then be realized by changing the 
"LO" phase d via the Larmor precession with compen- 
sation for optical shot noise as demonstrated by Juls- 
gaard et al. 20]. 

IV. RECONSTRUCTION OF 
SPIN-EXCITATION-NUMBER DISTRIBUTION 

The tomographic reconstruction of gives us an in- 
sight into the multi-particle entanglement of the spin 
systems when our two-mode (quantum and c-numbered 
"LO" modes) framework is employed. The concept of the 
multi-particle entanglement appears becasue of the fact 
that the basis states, e.g., |l)f (g) |7V— 1)^, are inherently 
entangled states as seen in Eq. Q. The complete infor- 
mation of density matrix can be obtained by the similar 
procedure of quantum optics outlined in Refs. [isL 12^12^ . 
The reconstruction procedure so simple as one restricted 
to just acquiring diagonal elements of the density matrix, 
which might be easier to implement, gives us clear phys- 
ical picture of multi-particle entanglement for collective 
spin states. Thus we now recapitulate a concrete proce- 
dure for determining the diagonal elements of the den- 
sity matrix in the entangled Dicke basis. By physically 
or numerically averaging the "LO" phase i? in Eq. lO, 
we have pj,|2l|23 

P(Q) = — / ddV{Q,d) - V 4^(Q) PMM, (12) 

where pmm = T (^^Ipt I-^)t ^^"^ 

(,) ^ ^ Af!(l-7?)^-^-"'r7"g^(g)exp(-Q^) 
* i'o {M-m)\m\ V^2"m! 

(13) 

with Hm{q) denoting the Hermite polynomials. Repeat- 
ing R times the measurement of Q affords the histogram 
{kv}, where represents the number of events that oc- 
curred in such a way that the measured value of Q be- 
longs to the bin with width 5Q and central point Q^. 
The probability of obtaining a specific histogram {k^} is 
given by 

p({fc.}i{PMAf})=i?!n^p^' (14) 

V 

where = P{Qy)SQ. Equations ((T^ and l(TI|l thus 
provide the relationship between the measurable his- 
togram {ky} and the distribution of spin excitation num- 
ber {pmm}- 

Two obstacles are encountered in reconstructing 
{pmm} from the measured histogram {k^}: the distri- 
bution {pmm} must satisfy the condition pmm > for 
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Spin-excitation Number 

(a) CSS 



Spin-excitation Number 



(b) SSS 



Spin-excitation Number 



(c) Dicke 



FIG. 1: (Color online) Simulated phase-averaged quadrature 
distributions {fci^}, where kv represents the number of events 
that occurred in such a way that the measured value of Q in 
Eq. belongs to the bin with width 5Q and central point Q^, 
(upper) and estimated spin-excitation-number distributions 
(lower, front row) for the (a) CSS |0>t, (b) SSS = 1)|0)t, 
and (c) Dicke state The true number distributions used 
to simulate the quadrature distributions {k^} are shown in the 
back row of the lower figures. The resultant cutoff parameters 
for minimizing the AIC (Eq. JTSJ) are (a) K = 1, (h) K = 9, 
and (c) K ^1. 



any M , and the number of unknown parameters {pmm} 
is large (i.e., the same as the atom number A'', although 
Pmm ~ for M » 1). To satisfy positivity constraints 
for {/Omm} and restrict the parameter space with a cutoff 
parameter a new set of parameters {9m} is defined by 



pii 



with e(^) = 6*^ + 6*^ H h 6*^. To search for the best 

estimate of {9m} and choose the most appropriate cut- 
off parameter K, the maximum likelihood method with 
Akaike's information criterion (AIC) is employed [25ll2^ . 
The best estimate of {9m} with the best cutoff parameter 
K can be found by numerically minimizing the AIC: 

A^''H{kA\{OM}) = -2 In [V^''H{K}\{ pmm})] + 2K. 

(15) 

Here, the first term denotes the log-likelihood function for 
the i^T-dimensional parametric model ofV {{k^}\{pM m}) , 
and the second term serves to suppress K and thus reduce 
the error due to redundant parameters [2^ |26| . 

Histograms {k^} for the CSS |0)t, SSS = l)|0)t 
(where S'(^) is the squeezing operator and Dicke 

state |l)f were generated by Monte Carlo simulation 
(20000 samples) with consideration of shot noise due to 
the optical probe. The histograms and reconstructed 
spin-excitation-number distributions for the three states 
are shown in Fig. ^ In this figure, the signal-to-noise 
ratio is fixed at 77/(1 — 77) = 1, since it can be esti- 
mated separately in experiments. The result for the SSS 
(Fig.n(b)) exhibits the collective-spin counterpart of the 



oscillatory number distribution The disappearance 
of odd-number spin excitations, such as and |3)-f, 
is a general feature of the squeezing below the standard 
quantum hmit 0. The phase-averaged SSS (Fig.n](b)) 
can thus be expressed by a classical mixture of a sepa- 
rable CSS and entangled Dicke states. In this way, the 
statistics of spin excitations in the entangled Dicke states 
give us clear physical picture of multi-particle entangle- 
ment for collective spin states [2^. It should be em- 
phasized that the characteristics of the small number of 
atoms ((a|a|) ~ 1) in the spin-up mode, i.e., the quan- 
tum mode, are emergent in the changes in the quadra- 
ture distributions via macroscopic number of atoms (e.g., 
(ajai) - 10" or lO^^ condensed in the spin- 
down mode, i.e., the "LO" mode. 



V. EXPERIMENTAL FEASIBILITY 

Recently the reconstruction of the CSS has been re- 
ported by Julsgaard et al. |^] in the context of the dis- 
persive quantum memory experiment. To observe the 
non-trivial spin-excitation- number distributions for e.g., 
the SSS and Dicke state, several schemes might be used: 

1) The aforementioned dispersive quantum memory 
scheme could lead to the observations of the spin- 
excitation-number distributions for the SSS and Dicke 
state by mapping a squeezed state of light and a single 
photon state, respectively, onto collective spin states. 

2) Instead of the coherent media transfer between light 
and coUactive atomic spin via dispersive quantum mem- 
ory, measurement-induced non-unitary schemes can also 
be utilized to generate the non-classical spin state. For 
the SSS the preparation and observation of the spin- 
excitation-number distribution shwon in Fig. ^ (b) is 
readily envisaged from the experimental results jlJl, l22| . 
To observe the Dicke state |1)|, the scheme proposed by 
Takahashi et al. ^2^ may be used by incorporating single 
photon detection. The so-called DLCZ scheme JOJ may 
also be used for generating the pseudo-spin analogue of 
the Dicke state, and the hyperfine-dependent scalar or 
tensor light-shift interaction Q can be used to mea- 
sure Q of Eq. © instead of vector light-shift interaction 
(Eq. dS))). Also noted is the experiment done by Black et 
al. 31] . With their cavity-based magnetic storage config- 
uration, the Dicke state generation and observation could 
be realized in a straightforward manner. 

3) Finally, we mention the resonant quantum memory 
scheme with dark-state polariton 32]. This scheme may 
also fit into our reconstruction scheme when one employs 
the experimental geometory similar to that of Black et 
al. 131|, i.e., a "coupling" laser beam and a "probe" laser 
beam for electromagnetically induced transparency form- 
ing the magnetic storage configuration. We think that it 
is worth doing to investigate further the difference be- 
tween the quantum states of collective atomic memony 
in dispersive scheme |20j and resonant scheme with dark- 
state polariton |3^ . 
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VI. CONCLUSION 

We show that the quantum state of the coUective spin 
excitation can be reconstructed from the quadrature dis- 
tribution, which can be obtained by the spin quantum 
nondemohtion measurements. The present reconstruc- 
tion scheme works for any states provided the mean spin 
excitation in the quantum mode is sufficiently small com- 
pared to that of the "LO" mode. Although this fact may 
exclude the application of our reconstruction scheme to 
e.g., the super-radiant Dicke state |iV/2)| (g) |A^/2)x [l3| 
and the cluster state {23|j the scheme may be a use- 
ful tool for investigating the multi-particle entanglement 
of atomic collective spin systems. Further extension 
to other systems such as collective nuclear spin excita- 
tion jsll can readily be envisioned. 
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N/2 



-N/2 



{tt/2) being the nontrivial part of the ro- 



taion matrix element for a spin-N/2 system 
the explicit expression [s^ ]. 



(-) 



m 



^(f +m)!(f -m)! ^ 



From 



(A3) 



and the de Moivre-Laplace theorem, we have a Gaussian 
distribution, given by 



Pz(m) 



exp[ 



N 
2 



(A4) 



as the large-N limit. This form is equivalent to the 
■y/ iV/2-foldcd quadrature distribution for the vacuum 
state, which is given by P^{q) =T('Zl(|0)TT(0|)k)T- 

As another example, the quadrature distribution for 
the Dicke state, = |1)tt(1| (|Jx = ~iN/2) + = 
— {N/2) + can be written as 



APPENDIX A: PHYSICAL MEANING OF SPIN 
QUADRATURE 

What is the physical meaning of the quadrature distri- 
bution for collective spin? Wc preliminarily set = 7r/2 
without loss of generality. The quadrature operator 
appears as a result of the approximation of in Eq. |0J. 
Thus the quadrature distribution P-\{q) (Eq. (fTUIl l is orig- 
inally 



N 

Pzim) = \{J,^m\J, = -- + !) 
and its large-N limit becomes 



Pz(m) 



(-2—) exph— ] 



(A5) 



(A6) 



Pz{m) = {Jz = toIpiIJ^ = m). 



Pz(to) = \{Jz^m\Jx 



(Al) 



For example, Pz(to) for the CSS, pt = |0)tT(0| {\X 
-N/2){Jx = -N/2\), can be written as 



(-) 

^2' 



(A2) 



This form is again equivalent to the a/ Af/2- folded 
quadrature distribution for the state of a single spin ex- 
citation, that is, P-\{q) = |(g|(|l)||(l|)|g)|. Thus the 
-y/iV/2-folded spin quadrature Pi{q) is just the proba- 
bility distribution of collective spin projecting onto the 
z-axis, i.e., Pz{fn), in the large-N limit. 



[1] W. M. Itano, J. C. Bergquist, J. J. Bollinger, J. M. Gilli- 

gan, D. J. Heinzen, F. L. Moore, M. G. Raizen, and 

D. J. Wineland, Phys. Rev. A 47, 3554 (1993). 
[2] F. T. Arecchi, E. Courtens, R. Gilmore, and H. Thomas, 

Phys. Rev. A 6, 2211 (1972). 
[3] M. Kitagawa and M. Ueda, Phys. Rev. A 47, 5138 (1993). 
[4] J. K. Korbicz, J. I. Cirac, and M. Lewenstein, 

Phys. Rev. Lett. 95, 120502 (2005). 
[5] V. Meyer, M. A. Rowe, D. Kielpinski, C. A. Sack- 

ett, W. M. Itano, C. Monroe, and D. J. Wineland, 

Phys. Rev. Lett. 86, 5870 (2001). 
[6] D. Oblak, P. G. Petrov, C. L. Garrido Alzar, W. Tittel, 

A. K. Vershovski, J. K. Mikkelsen, J. L. S0rensen, and 



E. S. Polzik, Phys. Rev. A 71, 043807 (2005). 
[7] S. Chaudhury, G. A. Smith, K. Schuiz, and P. S. Jessen, 

Phys. Rev. Lett. 96, 043001 (2006). 
[8] JM Geremia, J. K. Stockton, and H. Mabuchi, 

Phys. Rev. Lett. 94, 203002 (2005). 
[9] C. Orzel, A. K. Tuchman, M. L. Fenselau, M. Yasuda, 
and M. A. Kasevich, Science 291, 2386 (2001). 
[10] A. K. Tuchman, C. Orzel, A. Polkovnikov, and M. A. Ka- 
sevich, quant-ph/0504762 
[11] J. Hald, J. L. S0rensen, C. Schori, and E. S. Polzik, 

Phys. Rev. Lett. 83, 1319 (1999). 
[12] B. Julsgaard, A. Kozhekin, and E. S. Polzik, Na- 
ture (London) 413, 400 (2001). 



6 



[13] U. Leonhardt, Measuring the Quantum State of Light 
(Cambridge University Press, Cambridge, England, 
1997). 

[14] S. Schiller, G. Breitenbach, S. F. Pereira, T. Miiller, and 
J. Mlynek, Phys. Rev. Lett. 77, 2933 (1996). 

[15] A. I. Lvovsky, H. Hansen, T. Aichele, O. Benson, 
J. Mlynek, and S. Schiller, Phys. Rev. Lett. 87, 050402 
(2001). 

[16] B. Yurke, S. L. McCaU, and J. R. Klauder, 

Phys. Rev. A 33, 4033 (1986). 
[17] R. H. Dicke, Phys. Rev. 93, 99 (1954). 
[18] G. Klose, G. Smith, and P. S. Jessen, Phys. Rev. Lett. 86, 

4721 (2001). 

[19] A. Kuzmich, L. Mandel, J. Janis, Y. E. Young, R. Ejnis- 
man, and N. P. Bigelow, Phys. Rev. A 60, 2346 (1999). 

[20] B. Julsgaard, J. Sherson, J. I. Cirac, J. Fiurasek, and 
E. S. Polzik, Nature (London) 432, 482 (2004). 

[21] JM Geremia, J. K. Stockton, and H. Mabuchi, 
quant-ph/0501033 

[22] JM Geremia, J. K. Stockton, and H. Mabuchi, Sci- 
ence 304, 270 (2004). 

[23] K. Banaszek, Phys. Rev. A 57, 5013 (1998); K. Banaszek, 
G. M. D'Ariano, M. G. A. Paris, and M. F. Sacchi, 
Phys. Rev. A 61, 010304(R) (1999). 

[24] A. L Lvovsky, J. Opt. B: Quantum Semiclassical Opt. 6, 
S556 (2004). 

[25] H. Akaike, IEEE Trans. Autom. Control 19, 716 (1974). 
[26] K. Usami, Y. Nambu, Y. Tsuda, K. Matsumoto, and 
K. Nakamura, Phys. Rev. A 68, 022314 (2003). 



[27] W. Schleich and J. A. Wheeler, Nature (London) 326, 
574 (1987). 

[28] In Ref. the connection between the non-classicality 
of quantum state of light and the multi-particle entangle- 
ment of N qubit state was suggested. An observable crite- 
rion for non-classicality of light was proposed by W. Vo- 
gel, Phys. Rev. Lett. 84, 1849 (2000) and Th. Richter 
and W. Vogel, Phys. Rev. Lett. 89, 283601 (2002). Thus 
the sensitive entanglement testing may be done with this 
criterion. 

[29] Y. Takahashi, K. Honda, N. Tanaka, K. Toyoda, 
K. Ishikawa, and T. Yabuzaki, Phys. Rev. A 60, 4974 
(1999). 

[30] L. -M. Duan, M. D. Lukin, J. I. Cirac, and P. ZoUer, 

Nature (London) 414, 413 (2001). 
[31] A. T. Black, J. K. Thompson, and V. Vuletic, 

Phys. Rev. Lett. 95, 133601 (2005). 
[32] M. Fleischhauer and M. D. Lukin, Phys. Rev. Lett. 84, 

5094 (2000); Phys. Rev. A 65, 022314 (2002). 
[33] R. Raussendorf and H. J. Briegel, Phys. Rev. Lett. 86, 

5188 (2001). 

[34] R. N. Zare, Angular Momentum (John Wiley and Sons, 
New York, 1988). 

[35] J. M. Taylor, C. M. Marcus, and M. D. Lukin, 
Phys. Rev. Lett. 90, 206803 (2003); J. M. Taylor, 
A. Imamoglu, and M. D. Lukin, Phys. Rev. Lett. 91, 
246802 (2003); M. Eto, T. Ashiwa, and M. Murata, 
J. Phys. Soc. Jpn. 73, 307 (2004). 



